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We investigate the spinor solutions, the spectrum and the symmetry properties of a matrix- 
valued wave equation whose plane-wave solutions satisfy the superluminal (tachyonic) dispersion 
relation = ~ m? , where E is the energy, p is the spatial momentum, and m is the mass 
of the particle. The equation reads (i7^ dp, — 7* m)tp = 0, where 7* is the fifth current. The 
tachyonic equation is shown to be CV invariant, and T invariant. The tachyonic Hamiltonian 
= a ■ p + m breaks parity and is non-Hermitian but fulfills the pseudo-Hermitian property 
Hz{f) = P (—f) = V [rj'P^^ , where P is the parity matrix and P is the full parity 
transformation. The energy eigenvalues and eigenvectors describe a continuous spectrum of plane- 
wave solutions (which correspond to real eigenvalues for |p| > m) and evanescent waves, which 
constitute resonances and antiresonances with complex-conjugate pairs of resonance eigenvalues 
(for IpI < m) . In view of additional algebraic properties of the Hamiltonian which supplement the 
pseudo-Hermiticity, the existence of a resonance energy eigenvalues E implies that E* , —E, and 
—E* also constitute resonance energies of H5. 
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I. INTRODUCTION 



Extensions of the Dirac equation have been discussed in great numbers in the scientific hterature, including field- 
theoretical formulations [1] . Notably, there have been efforts to incorporate space-time curvature into the equation [2] , 
and, to enhance the equation for spin-1 particles such as the photon, and to map the Maxwell equations onto Dirac- 
type equations [3]. Here, we follow a different route and analyze whether it is meaningful and consistent to modify the 
mass term in the free Dirac equation so that the energy-momentum dispersion relation changes from E'^ = p ^ + 
to E'^ ~ p"^ — . The latter dispersion relation would be relevant for superluminal particles, which travel faster than 
the speed of light. The tachyonic spin-1/2 equation reads 

(i7^ - 7' rn) xl>{x) = , (1) 

where the 7'^ arc the Dirac matrices, = djdx^ is the derivative with respect to space- time coordinates, and 
7^ = 17° 7"'^ 7^ 7*^ is the fifth current [here, x — {t,f) denotes a space-time four vector]. Because of the occurrence of 
the Dirac matrices, we propose the name 'tachyonic Dirac equation'. The corresponding Lagrangian density reads 

£ = ^V^7^ - V = ^ (^7' 7^(5^ V-) - 7' 7^^) - mi^i; , (2) 

where = ip+ 7" is the Dirac adjoint. It is perhaps not immediately obvious that the 7^ matrix has to be placed 
with the current matrix 7^ in the Lagrangian; otherwise an inconsistent equation is obtained for ip upon variation 
with respect to ^. 

Indeed, the tachyonic Dirac equation has been indicated by Chodos, Hauser and Kostelecky in Refs. [4] in the 
equivalent form (17^ 7'^ df^—m)tp = 0, where the current is modified instead of the mass term (see also Refs. [5, 6]). More 
recently, the noncovariant (Hamiltonian) form of the equation has been indicated in Refs. [7, 8], with a preliminary 
ansatz for its solution. Very recently, the equation has been used in the context of the Gross- Neveu model in Ref. [9] . 
However, neither the complete set of solutions, nor the spectrum, nor the symmetries of the equation have been fully 
investigated in the literature to the best of our knowledge. 

An important observation made in the current work concerns the Hamiltonian 

H5 = a-p + l3-f^m (3) 

which is obtained from the tachyonic Dirac equation. Here, a = 7*^7 and /3 = 7*^, with the identification of the 
Dirac matrices discussed below in Eq. (6). The Hamiltonian breaks parity and is non-Hermitian but fulfills the 
pseudo-Hermitian equation 



(4) 
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where P is the parity matrix and V is the full parity transformation. While the scalar product induced by V is not 
positive semi-definite, a few essential properties of the Hamiltonian persist under the pseudo-Hermiticity and help to 
characterize the spectrum. We thus extend recent work on non-Hermitian, but pseudo-Hermitian or PT-symmetric 
Hamiltonians [10-20] which have also been applied in the field-theoretical context [21-23]. 

We proceed as follows: In Sec. II, we analyze the transformation properties of the covariant formulation of the 
tachyonic Dirac equation under charge conjugation, parity transformation and time reversal. The propagator and the 
explicit solutions of the tachyonic Dirac equation arc analyzed in Sec. III. Finally, conclusions arc drawn in Sec. IV. 
Units with h = c = cq = 1 are used throughout the paper. 

II. DISCRETE SYMMETRIES: C, V, AND T 
A. Charge Conjugation 

We start from 

[r{\d,,~eA^)~j''m]^ = 0, (5) 
where the are the Dirac 7 matrices, which we use in the Dirac representation. 

The Dirac 7*^ matrices fulfill the following anticommutator relations, 

{7^7'^} = 2g^^ (7) 
where g^'^ = diag(l, —1, —1, —1) is the space-time metric. Transposition and complex conjugation leaves 7^ invariant, 

V'+((7'')+(-iK-eA^) =0, (8) 

where the partial derivative acts on the left. The Dirac adjoint is introduced as 

(^+7°) 7° ((7")+ (-i^M - e A,,) - 7^™) 7" = . (9) 
Using the identities 7° (7^)''' 7° = 7^ and ■y^j^j^ — —7^, one writes 

V;(7^(-iK-e^p) +7'™) =0. (10) 

A further transposition leads to 

((7^)^(-ia^-eA^)+75m) =0. (11) 
Introducing the charge conjugation matrix C with the property 

cirfc-'^-r, (12) 

the charge conjugation of the fifth current is calculated symbolically as 

C (7^)^ = Ci{j'f {j'f (7^)^ (7°)^ 

= iC (7')^ c-1 C (7')^ c (7')^ C (7°)^ 

= i (-7^) (-7^) (-7^) (-7°) 

= 17^72717" = 1707172^3 =75. (13) 

Therefore, 4'^ = C ip^ fulfills the charge conjugate tachyonic Dirac equation, 

[7^ (i + e A^) + 7^771] ^.'^ = . (14) 

As a result of the calculation of the Dirac adjoint (7^ ^"--yS-yO _ —7'^) and the charge conjugation using the C 
matrix, the term with the 7^ matrix has changed sign. In contrast to the ordinary Dirac equation [24] , the tachyonic 
equation is not C invariant. 



3 



B. Parity 

Again, we start from 

(ij^d^, - 7'5m) i^{x) = (i7°9o + - = . (15) 

Under parity, x ^ Xp = {t, — f) and so 

(i7"ao-if 5,-7'™) V-C^p) =0. (16) 
One introduces a parity transformation matrix P, where P = 7° in the Dirac representation, with the property 

Pj°p-^ = 7" , Pj'p-^ = -Y , Pj^p-^ = -7^ , (17) 
where again i = 1, 2, 3 is spatial. So, ijp^ [x) = Pip{xp) fulfiUs the parity transformed tachyonic Dirac equation, 

P {i-f"dQ - i7*9, - 7^ m) ^/.^(.t) = (18) 

which can be rewritten as 

(i7°ao + il'd, + 7^ to) (x) = . (19) 

In covariant notation, we have 

(i7^ + 7^m) (x) = . (20) 
The tachyonic Dirac equation is not parity invariant. 

C. Time Reversal 

Starting from the tachyonic Dirac equation, 

(i7^c>M-7'm) ^ = 0, (21) 

transposition and complex conjugation leads to 

^+(-i(7^)+K-7'm) =0, (22) 

where again the differential operator d ^ acts to the left. The insertion of the 7'^ matrix ensures Lorentz covariance 
and leads to 

(^+7°) 7O (-i (7^)+ - 75 to) 7° = , (23) 

which is equivalent to 

rp (-i7^ + 7'^ TO.) = . (24) 

We here use the fact that a proper formulation of time reversal must respect the fact that the in and out states 
interchange under T. We thus have to formulate the time reversal in terms of the adjoint of the Dirac equation (see 
Sec. 4.4.2 of Ref. [25]). Again, we transpose and write 

i (7^)"^ + 7^ to) i/i^ = . (25) 

We introduce the time reversal matrix T, with the properties (T means time reversal, T is the time reversal matrix, 
and T denotes the transpose), 

T(7")^T-i =7", r(f)^T-^ = -f. (26) 
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The time reversal of the fifth current is calculated as 



= iT (7^)^ T (7')^ T-i T (7')^ T (7")^ 



= i(-7^) (-7^) (-7^) (+7") 

= -17^727^° - -i7°7S^7^ = -7^ • (27) 
The time-reversed tachyonic Dirac equation thus reads as 

(i7^ - 75 m) V'^(.t) = , (28) 

where iji''{x) = T ^p^ {xt) with x = {—t,r). As a result of the replacements 7'^ — s- 7°^^^" and 7^ ^ Tj^T^^, the 
tachyonic Dirac equation is seen to be T invariant. 

D. Overview of the Symmetry Properties 

A brief summary is in order. The charge conjugation matrix fulfills [24, 26] 

C(7^)'^C-i = -7^. (29) 
In the Dirac representation, a possible choice is C = 17^ 7". The parity transformation fulfills 

P^op-i ^ ^0 ^ pyp-i = _y . (30) 

In the Dirac representation, we may choose P = 7". The time reversal operation fulfills 

T(7")'^T-i =7", T{YyT-^^^Y. (31) 

In the Dirac representation of the 7 matrices, one can choose T = 17^ 7^. The matrices C, P and T fulfill the relations 
= C^^, = P^^, = T^^, where by T we denote the transpose. A representation of CVT is given by 

CPT = 7^ . (32) 

This implies that 

{CPT)-f''{CPT)-^ = --fi" , (33) 

a relation which has been given on page 239 of the textbook [27] and in Eq. (40.48) of Ref. [28]. We also recall the 
transformation properties of 7^ just derived, 

C(7^)^C7-l=7^ 

P-f^p-^ = -7^ , (34) 

^(7^)^^-l = -7^ 

We note that the transformation property for T depends on our precise formulation of the time reversal operation as 
detailed in Sec II C. According to Eqs. (14), (20), and (28), the tachyonic Dirac equation is separately CV invariant, 
and T invariant. 

III. SOLUTIONS OF THE TACHYONIC EQUATION 
A. Plane— Wave Solutions and Chirality 

We start from the eigenfunctions of the operator a ■ k, which are given by 



'(f) 



cos TT \ / — sm c 



-lip 



".+(^=1 .7;p.i. -wH i. (^^a) 
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and fulfill 

(7 • k 



\k\ 



a±{k) = ±a±{k) , (35b) 



Va±(fc)®ai(fc) = l2, y2i-lf a±{k) ® a+{k) ^ ^ , (35c) 
± ± 1^1 

where 9 and ip are the polar and azimuthal angles of the wave vector k. The normalized positive-energy chirality and 
helicity eigenstates of the massless Dirac equation are 

These eigenstates immediately lead to plane- wave solutions of the massless Dirac equation, which are also eigenstates 
of the chirality 7^ and of the helicity E • p/\p\, where E is the vector of (4 x 4)-spin matrices, 

V'(a;) =M±(fc)e-"'-^, k = {E,k) , E = \k\ , (37a) 
i-,^ d^iPix) = k^,pix) = , ido^Pix) = E,pix), (37b) 

^u^{k) = ^^u^{k)^±u^(k), ^«±(fc)c'^'- = ±j.±(fc)e'^^-. (37c) 

\p\ 

For the negative-energy solutions of the massless Dirac equation, we need 

„,m.J^ ._(?, V (38) 

The negative-energy states fulfill the relations 

0(x) = v±{k) e"'"^ , k^{E,k), E ^ \k\ , (39a) 
i 7^ a^0(x) = - 7^ k^^{x) = , i do<j)ix) = -E 0(x) , (39b) 

^ v±{k) = v±{k) ^ ±v±{k) , ^^;±(fc)e-'^'- = T^;±(fc)e-'^^•^ (39c) 

The subscript of the v vectors is chosen according to the corresponding eigenvalue of the chirality 7^, i.e., 7^ u± = ±m±, 
and 7^ w± = ±z;±. For negative-energy solutions, chirality and helicity are opposite [24]. It is worth noting that the 
set of four states, 

u±(A?)e*^"', v±i~k)c'^-^, (40) 
form a complete and orthonormal set of eigenstates of the massless Hamiltonian Hq, which is a Hcrmitian operator, 

Ho = d-p, a=:7°7, (41a) 

Fou±(fc)e'^-*'= \k\u±{k)e'''-^ , Ho w± (-fc) c'^'*" = -|fc| w±(-fc) e'^■^^ (41b) 

Ua{k)ufi{k) = 6ai3 , v+{-k)v/3{-k) ^ Safs , u+{k)vfj{-k) ^ , (41c) 

where a, f3 = ±. These relations hold in view of a'\_{k) a^{-~k) = al(A:)a_(— fc) = 0. In view of the relation 
{j^kfi — 7^™)^ = k"^ + m? = E'^ ~ k'^ + m? , the generalization of these solutions to the massive tachyonic Dirac 
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equation is rather straightforward. We find 

u+ik) -~ 



7^ TO — 7^ 



V2 ^ (E ~ \k\)^ + 
( m-E+\k\ 



■.u+ik) 



V2^jiE-\k\rj 

m + E - \k\ 



for positive chirahty and 



U-ik) 



\ V2^iE-\k\r 



7^ k^ — m 



a+{k) 



a+ik) 



(42a) 



^/2^(^;-|fc|)2 + m2 
/ m + E-\k\ 



V2 (E - \k\)^ + 
-m + E- \k\ 



a-{k) 



\ V2^JiE~\k\)^+m^ 



a_(fc) 



(42b) 



J 



for negative chirality. Here, E = v^fc^ — m? with k^ > m^. For to — > 0, the particle has to be on the mass sheh, 
therefore E — > |fc|, and we have the tachyon spinors approaching their massless limit, i.e., U+{k) u+(fc) and 
U-{k) — > u-{k). The negative-energy eigenstates are 



V+{k) 



-f^m + jf^k^ 
, u+{k) 

V2JiE-\k\)^+m^ 



( -m-E+\k\ 



^J [E - \k\f + m? 
-TO + E -\k\ 

\^ V2\/{E- |fc|)2 +to2 



a+{k) 
a+{k) 



(43a) 



for positive chirality (negative helicity) and 

V-ik) = 



—7^ fc^. — 7^ m 



V2 ^ (E ~ \k\)^ + 
( -m + E-\k\ 



^^{E-\k\f+: 
m + E - \k\ 

\ V2^iE-\k\y+: 



-ik) 



a^{k) 



a^{k) 



(43b) 



for negative chirality (positive helicity). For m — > and E — > |fc|, we have V+(fc) — > w+(fc) and VL(fc) — > v_(fc). These 
states are normalized, i.e., UX{k)U+{k) = Ut{k)U-{k) = VX{k)V+{k) ^ V^{k)V-{k) = 1. The corresponding 
positive-energy solutions of the massive tachyonic Dirac equation arc given as 



*(a;) = C/±(fc)e-"=-^, k^{E,k), E 



fc2 



(i 7^ df, - 7 V) «'(x) = (7^ /.V - 7^71)1' (x) = , 



(44a) 
(44b) 
(44c) 
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FIG. 1: (Color.) Illustration of the velocity addition theorem f{u,v) = {u + v)/{l + uv), in the su- 
perluminal domain with it G ( — 1, 1) and v £ (1, 3). The superluminal velocity v remains superluminal 
upon transformation from a frame which moves at subluminal velocity u with respect to the observer; 
the domain —1 < f(u,v) < 1 is excluded, consistent with Eq. (63). 



The negative-energy solutions are given by 

^{x) = V±{k)e'^-'' , k^{E,k), E=\lp~m'^ , (45a) 

(17^ - 7^m) $(x) = (-7^' - 7^n)$(a;) = , (45b) 

\da<^{x) ^ -E<i>{x) , (45c) 

All of the discussed solutions fulfill the superluminal dispersion relation 

= P - (46) 

and are eigenstates of the tachyonic Dirac Hamiltonian (49), as given below. 

B. Plane— Wave Solutions and Hamiltonian 

It is instructive to investigate the Hamiltonian form of the tachyonic Dirac equation, which reads 

i7^^-7'"i^ ^(.t) =0. (47) 

Multiplication by 7*^ leads to the non-covariant Hamiltonian formulation, 

H^iPir)^ {a-p + P"f^m) ij{T^= E^Pir), (48) 
where c5 = 7^^ 7 and = 7*^. The tachyonic Dirac Hamiltonian obviously reads as 

Hz = H^if) = a-p + /3j^m, (49) 
where p = — iV is the momentum operator. The Hermitian adjoint is 

H+{r)^d-p-Pj^m, (50) 
where by the superscript + we denote the Hermitian adjungation. Therefore, 

PH+{f)p-^ = P {a ■ p - p -i^ m) P-^ 

— — d ■ p + m , (51) 
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because (f^j^)'^ = (7^)''' (7°)''' =7^7° = —7^7^. Here, p —>■ = p because the Hcrmitian adjoint of the 
momentuin operator is the momentum operator itself. It is thus easy to check that the tachyonic Hamihonian Hc^ is 
a pseudo-Hermitian operator [15-18, 29], with the property 

i75(r) = PH+{-r)p-' = V H+{r)V-^ , (52) 

where again P = 7*^ is the matrix representation of parity and V is the fuh parity transformation. Hamiltonians with 
this property have been analyzed about seventy years ago (see the work of Pauli, Ref. [29]), and have been the subject 
of recent investigations [10-12, 14-18] because under rather general assumptions, they have real energy eigenvalues, 
or resonance eigenvalues which come in complex-conjugate pairs. 

This can be seen as follows and holds even though the the scalar produce (V'l V'2) is not positive-semidefinite. First 
we recall that because the spectrum of a Hermitian adjoint operator consists of the complex conjugate eigenvalues, 
there will be an eigenvector (pif) with eigenvalue E* provided there exists an eigenvector il){r) with eigenvalue i?, 

i/g (f ) 7/. (f ) = P 7/. (f) , H+ (f) 0(f) ^ E* 0(f) . (53) 

The transformation f — > —r and the insertion of the parity matrix P = 7° leads to 

H+{-r)c^{-r)^E*<P{-r), 

PH+{-r)p-^ {P(j){-f)) = E* P(j){-f) . (54) 

By assumption, PH^ {—f^P"^ = Hc,{r) and thus 

H^(f)P^(-r)^E* Pcj)i-r), 

Hsir) ipir) = E*ip{r) , ip{f) ^ P<f){-r) . (55) 

This implies that ?/'(f) = P(f)[-~r) is an eigenvector with eigenvalue E* . The eigenvalues of i/s thus come in complex- 
conjugate pairs. Furthermore, an additional quasi-pseudo-Hermitian property can be established, 

H^{f) = -pH+{f)p-\ p = j"^', P=-p-'. (56) 

Going through the above derivation again, we easily find that if E and E* arc energy eigenvalues, then —E and 
—E* will also be energy eigenvalues. This configuration of eigenvalues can be characterized as a "St. Andrew's cross 
configuration" in the complex plane (see also Ref. [30]). 

For a number of operators such as the VT symmetric variants of odd anharmonic oscillators [13, 19, 20], the 
eigenvalues are all real. However, for the tachyonic Dirac Hamiltonian, we shall need the full spectrum which consists 
of both real as well as complex resonance eigenvalues. In order to characterize the spectrum of resonances, it is 
instructive to show that the resonance energies of the tachyonic Dirac Hamiltonian are real in the second-order 
formalism. A solution of the tachyonic equation 

i7^-^-7^m) V(a;) =0, ^/'(a;) = f/ e""^-^ , (57) 



with k = (E, k), must also fulfill the Hamiltonian eigenvalue equation 

(d-k + Pj^rnj U = EU , (58) 

and for "ipi^) = V c^*''^ we have 

[-d-k + I3^^v^ V = -EV . (59) 

Squaring these equations, we thus have 

(a • fc -f /37^ m)^ C/ = (fc^ _ „^2^ U = E^ U , (60) 

and 

{-a-k + p-i^m^^ V ={k'^ -rri^^V = E'^V . (61) 
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We here observe that the matrix /37^ is a representation of the imaginary unit, in view of the identity {P^^Y = — 14x4. 
Connections to other representations of the imaginary unit within superluminal spin-1/2 Hamihonians are discussed 
in A. The superluminal energy-momentum dispersion relation 

- (62) 

(with real E) is thus fulfilled. Resonance energies E of the tachyonic Dirac Hamiltonian thus have to be either 
completely real (vanishing imaginary part) or completely imaginary (vanishing real part), so that E'^ is real. 



Need for the Inclusion of Resonances and Anti— Resonances 



One might wonder about the nature of plane-wave solutions with |fc| < to, which according to E = — 
would have imaginary energies. Denoting by it > 1 the velocity of the particle (expressed in units of the speed of light 
c), we have for tachyonic particles, according to Refs. [31-35], the relations 



E = 



1 



|fc| = b1 = 



1 



> m . 



(63) 



where u is the particle velocity in units of the velocity of light, c. Superluminal Lorentz transformations have been 
discussed in Refs. [31, 36] (see also Fig. 1); it is impossible to 'stop' a superluminal particle, and the 'rest frame' of a 
superluminal particle is that of infinite velocity [32, 37]. We conclude that since |fc| > to, the problematic imaginary 
energies which are encountered for |fc| < to are excluded by the tachyonic dispersion relation. This leads to a restriction 
|fc| > TO for the permissible values in constructing stable tachyonic wave packets ("Fourier transformations") from the 
fundamental plane-wave solutions. In a distant analogy to the "positive definite" functions which are the subject of 
Bochner's theorem (Theorem IX. 9 of Ref. [38]), the wave packets of tachyonic particles have to be "positive definite" 
in the sense that they are naturally constructed using plane waves with |fc| > to in view of the tachyonic dispersion 
relation. (A "positive definite" function has a positive Fourier transform, or more precisely, is characterized by a 
Fourier transform which constitutes a positive Borel measure.) 

In Refs. [31-34], considerable effort has thus been invested into the elimination of the states with < m? from 
the scalar tachyonic theory, because it was realized that the resonances correspond to evanescent waves. However, we 
here find that these states qualify themselves as resonances of the tachyonic Dirac Hamiltonian. 

We define the width 7 of a resonance of the tachyonic Dirac Hamiltonian as follows, 



= ± \/fc2 - ,7i2 -ie = =Fi- , r = 2YTO2-A;2^ 
The wave functions describing the resonances are as follows. 



(64) 



R-{k) 





TO -I- 


+ 
2^ ^ 


\k\ 








-|- to2 


+ \ 






TO — 


2^ 


\k\ 




a/2 




-|- to2 


+ 1 


Y2 




TO — 


2^ 


\k\ 




V2 


7^2 


+ to2 


+ \ 


Y2 




-TO - 


2^ 


-\k\ 




V2 


Jk^ 


+ to2 


+ \ 


Y2 



a+{k) 



\ 



a+{k) 



a^{k) 



a^{k) 



(65a) 



(65b) 



E^ - iF = -iy^TO-P, P<m^. (65c) 
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The antircsonance eigenstates are 



S+{k) 



S-{k) 



E 



/ - 


-m — 


2^ ^ 


-\k\ 






j 

JP 




+ i 






-m + 


2^ 


-lfcl 












Y2 


( - 


-m + 


2^ 












1 

+ 4 


p2 


m + 


2^ 


\k\ 






A^ 


+ 


+ 1 




ir = ^ 

2 2 


\J m 


- fc2 


1 


k 



a+{k) 
a+{k) 

a^{k) 
a-{k) 



\ 



(66a) 



k < m 



(66b) 



(66c) 



These states are normalized, i.e., R+{k)R+{k) = Rt{k)R-{k) = S'+(fc)S'+(fc) = S^{k) S-{k) = 1. For |fc| 
have E' — > 0, and the eigenstates read as 



U+{k) R+{k) - 

U_{k) R_{k) ^ 
The negative-energy spinors tend to the following values, 

v+{k) ^ s+{k) ^ 1^ 

V-{k) S-{k) { ' 



a+{k) \ 

/ 


-a_(fc) 



\k\ 





-a+{k) 

-a^{k) 




m , 



(67a) 
(67b) 

(68a) 



\k\ 



According to Eq. (63), these spinors, which correspond to a tachyon of infinite speed, correspond to the solutions 
which describe a tardyon spinor at rest. This point has also been made in Ref. [37]. 
At time t = 0, the above vectors are eigenstates of the Hamiltonian, 

(69a) 
(69b) 
(69c) 
(69d) 

with k E ]R^. These form a complete set of eigenstates of the tachyonic Dirac Hamiltonian and allow us to solve the 
initial value problem. 





= C/±(^)e'^'•^ 


|fc| > m , 




= 14(-fc)e'''•^ 


|fc| > m 




= i?±(fc)e'^'•^ 


fc < m , 




= 5±(-fc)e'^'•^ 


< m 



idtij{t,r) = H5^{t,r), 7A(0,r)-*(r=), 



(70) 



where the initial wave- function configuration is given by ^(r). We should indicate a little caveat, here: The U±{k) 
are mutually orthogonal (positive energy), even for the same fc, and the V±{k) (solutions for negative energy) have 
the same property [see Eq. (76c) below]. The initial value problem is thus solved on the level of the relativistic 
quantum theory for both positive-energy as well as negative-energy wave packets. The only remaining problem is the 
non-unitary character of the time evolution generated by the tachyonic Dirac Hamiltonian through the resonances, 
and antiresonances. Such a non-unitary time-evolution in the presence of resonances is discussed in detail in Ref. [39] 
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for the case of the odd anharmonic oscihator. We should clarify that the resonance solutions R±{k) and S±{k) are not 
resonances in the sense of a virtual state in a cross section (which "decays" into a final state) , but constitute complex- 
valued resonance energies of the Hamiltonian matrix. They describe the evanescent subluminal components of a wave 
packet which is genuinely superluminal. According to the ie prescription from Eq. (64), resonance solutions constitute 
the analytic continuation of positive-energy solutions for |fc| < m. They are exponentially damped for propagation into 
the future. By contrast, anti-resonance solutions describe the continuation of negative-energy solutions for |fc| < m, 
these are damped for propagation into the past. This prescription ensures that an inconsistent exponential growth of 
the wave packets cannot occur [1] . 

Our complete set of functions (not restricted to |fc| < m) allows us to represent localized wave packets in terms of 
a Fourier decomposition. One example of a localized wave function in a scalar theory is 

It has the property 

' d^r'S^{r')S^{f+r')=5^{r}. (72) 



Likewise the function 



has the property 



7:ri3^+(fcK''^+ /- Tr^U+{k)e''^-^ (73) 

J\k\<7n (27r)-^ J\k\>rn [^■^) 



^ ^h.i^ f di k 



dV^+{^-^)^{r + r')= J ^ ^{k) e"'-^= J ^ c"^-~ = 6\r) . (74) 

where ^(fc) is the Fourier transform of 'i>{f). 

Three remarks are in order, (i) In an idealized position measurement using a detector, according to the wave 
function collapse postulate of quantum theory, the wave function is proportional to a Dirac-5 after the measurement. 
It is important to be able to localize a tachyonic particle in order to measure it. According to the Heisenberg principle, 
the position of quantum particle cannot be measured significantly better than its Compton wavelength, but this finding 
does not imply that idealized position measurements are meaningless. In our formalism (see also Ref. [1]), we prefer to 
retain localizability of a particle rather than prevent the emergence of the small region of non-unitary time evolution in 
the sector |fc| < m. (ii) It would be somewhat unnatural to exclude the resonances and antiresonances of the tachyonic 
Dirac Hamiltonian from the discussion a fortiori. Canonically, resonances and antiresonances play an important role 
in the spectrum of quantum Hamiltonians and should receive a physical interpretation [40]. In our case, we propose 
an interpretation in terms of the evanescent waves that describe tardyonic components of a genuinely tachyonic wave 
packet. ( Hi) For completeness and for general interest, we note that the Hamiltonian H' = a- p + i /? m, where i is the 
imaginary unit and m is a real mass term, is connected with by a unitary transformation, which is constructed 
in A. 

With regard to remark (i) and Refs. [31-35] , we observe that full unitarity cannot be preserved anyway in a tachyonic 
theory at the one-loop level, as observed in Ref. [41]. Although the unitarity violation at the one-loop level could 
be interpreted as an argument against the consistency of tachyonic field theories, we here prefer to countenance a 
small unitary violation [1] due to the evanescent waves than to sacrifice the localizability of a wave packet. Likewise, 
in Ref. [42], it was argued that tachyonic particles cannot be described by a unitary representation of the Lorentz 
group. In general, we have to sacrifice either unitarity or localizability in we wish to incorporate tachyons into field 
theory. As observed above and explained in Ref. [1], the slight violation of unitarity is restricted to a relatively small 
kinematic region. Upon the inclusion of the resonances, we observe that one can write field commutators like those 
encountered in Eq. (1.4) of Ref. [34] in terms of the full Dirac-5 function instead of the "filtered" Dirac-(5 function 
given in Eq. (1.5) of Ref. [34], which is not localized. This is explained in detail in Ref. [1]. 

With regard to remark (iii), we note that the field theory based on the Hamiltonian H' was analyzed in Ref. [37]; 
this analysis has recently been enhanced in Ref. [30]. 



D. Orthogonality Properties of the Solutions 



Finally, a discussion of the orthogonality properties of the solutions of the tachyonic Dirac equation is in order. 
Let us suppose that 'ip2 is an eigenvector of a pscudo-Hermitian Hamiltonian H5 with eigenvalues E2 and tpi is an 
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eigenvector of with eigenvalues Ei . Then, 

E2 {'^l\V\^2) = {^1\VE2\'<}J2) = (V'l|Pi?5|V'2) 

= {^i\VH^V-^V\iJ2) = {i^i\H+V\^P2) 

= (i?5V'i|7'|^2> = El (Vji|7'|V2> , (75) 

i.e., eigenvectors for different energy are orthogonal with respect to the T'-scalar product. Our plane-wave solution 
form a complete set of eigenstates of the Hamiltonian i/s, with the properties 

iJs U± (k) e'''-^ = EU± (k) (^^-^ , (76a) 

V± i-k) e'*^-"^ = ~EV±{~k) e'^-"^ , (76b) 

U+ik) Up{k) = , V+{-k) Vp{-k) = S^fi , (76c) 

C^^(fc)7°^/3(+fc) =0, a,/3 = ±, (76d) 



and E = \/k'^ — ra?- . The relation J7+(fc)7''Va(+fc) = is equivalent to 

{u±(k)(^^-^\V\ y±(fc)c-'*^""^ = l^±i$)^^'"\f\ Y±{^k)(f-^'^ = 0, (77) 

which is equivalent to the orthogonality of the eigenvectors with different energies E and —E under the scalar 
product defined in Eq. (75). Within the space of positive-energy and negative-energy solutions, the time evolution 
is unitary with respect to the ordinary scalar product defined in Eq. (76c). Furthermore, we observe that while in 
general C/+(A:)V/3(— fc) 7^ 0, the four solutions U±{k) e''^ '" and V±(— fc) e''^''" form a linearly independent set for given k. 
However, as explained in B, the V norm of the plane-wave states with a nonvanishing k actually vanishes, and the V 
scalar product therefore is not positive definite. Therefore, as explained in detail in Ref. [43], it would be desirable 
to work with a different scalar product, obtained by prepending the V operator with a second operator, called C 
in Ref. [43] (while being manifestly different from charge conjugation, see B), so that the ensuing redefined scalar 
product is positive-definite. 

The construction of a positive-definite scalar product, conserved under the pseudo-Herniitian time evolution, is 
left as an open problem for future investigations. However, one important observation can be made at the current 
stage: namely, the time evolution induced by the pseudo-Hermitian Hamiltonian is separately unitary for wave packets 
constructed from exclusively positive-energy and negative-energy wave packets, in view of the property (76c) and the 
reality of the eigenvalues. This pseudo-unitarity is sufficient if one assumes that the Hamiltonian naturally breaks up 
into a positive-energy and negative-energy Hamiltonian, according to the interpretation of positive-energy solutions 
describing particles and negative-energy solutions describing antiparticles [1]. Indeed, for a Hermitian Hamiltonian, 
the scalar product 

{Mt)\Mt)) = J d''xi'+ir,t)^{r,t) (78a) 

is conserved, as can be seen by differentiation with respect to time. For a pseudo-Hermitian Hamiltonian which fulfills 
Eq. (52), time evolution is unitary with respect to the scalar product (i/'i(t)|7'|V'2(i)) (see Ref. [29]). This is because 



[dtiMtmMt]) = J d'.^ (-iatVi(i,^)+ PV'2(i,r) + J d^xi4{t,r)VidtMt,r) 
d-\ (-i75(r)Vi(t,r))+ rMt,r)+ f d^x V^+(t, r) i?5(r)^2(i, r) 



{Mt)\H+irm2) + ii'ilV H5(r)|V^2) = , (78b) 
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where we use the property H^{r)V = V H5{r). Thus, (?/'i(t)|'P|'02(i)) is conserved under pseudo-Hermitian time evo- 
lution. However, for wave packets constructed from plane waves that fulfill the superluminal dispersion relation (63), 

we can actually show a stronger statement because the energy eigenvalues Et = — m? with > m are purely 
real. Let us study two positive-energy normalized wave packets of the form {j ~ 1,2), 

V',(i,r) = j (w,Ak)U+{k)+w,Ak)U-{^)) c-^^^'+'^^, 

J d^k {\w,,+ ik)\' + h,,-(fc)p) = 1 , (78c) 
where the U±{k) are given in Eq. (42). Then, 

d'^y y {wl4k')Ut{k') + wl_{k')U+{k')) 

X {w2.+ ik")U+ik") + W2Ak")U-{k")) c-"^Ez"-Ez,)t+Kk"-k'yr 
= /d'fc" {wl+{k')w2Ak")Ut{k')U+{k") 

+wl _{k')w2.-{k") Ul{k') U-{k")^ c^'(-^E"-^£''*(53(fc" - k') 

^ j K+(^)^2,+(fc)+<_(fc)^2.-(fc)) 

= (Vi(O,r)|7^2(O,?0)- (78d) 

Thus, unitary time evolution holds provided the energies — fulfill the dispersion relation given in 

Eq. (63), which implies that they are manifestly real. Within the manifolds of positive-energy and negative-energy 
wave packets, the C operator in the sense of B thus can be taken to be equal to the parity operator V, in accordance 
with the findings of Ref. [44] . 



IV. CONCLUSIONS 



In the current article, the tachyonic Dirac equation (17'' 5^ — 7^771) -0 = is analyzed. We show that it is invariant 
under time reversal (T, Sec. II C), and invariant under charge conjugation and parity (CV). Our formulation of the 
time reversal operation involves the Dirac adjoint (see Sec. 4.4.2 of Ref. [25]). We generalize the treatment of the C, 
V and T transformations to the tachyonic Dirac equation in Sees. II A, II B and II C. These are the desired symmetry 
properties for the description of neutrinos. The solutions of the tachyonic equation approximate the chiral solutions 
of the massless Dirac equation in the limit m — >■ (see Sec. Ill A). In writing the solutions, one is naturally led to 
the use of a basis of chiral eigenfunctions, aligned with the propagation four-momentum k = {E, k) of the particle. 
The tachyonic Hamiltonian given in Eq. (49), = a ■ p + P m, is shown to be pseudo-Hermitian. The concept 
of pseudo-Hcrmiticity is related to an article of Pauli [29] who established Hamiltonians with the given property 
as viable alternatives to Hermitian Hamiltonians (see also Refs. [10-12, 14]). In particular, pseudo-Hermiticity has 
been indispensable in the recent generalization of the so-called Bender-Wu formulas [45-47] to odd anharmonic 
osciUators [19, 20]. 

In the tachyonic Dirac Hamiltonian iJs = a ■ p + /Sj^m, the matrix 14 = /3 7^ is a representation of the imaginary 
unit in four dimensions, 

where l2x2 denotes the 2x2 unit matrix, and lixi denotes the 4x4 imit matrix. Indeed, it was Dirac who 
argued [48, 49] that the linearization of the Klein-Gordon equation required the introduction of 4 x 4 matrices; the 
most straightforward implementation of the imaginary unit in four-dimensional space is due to the matrix i4 = /? 7^ 
because it generalizes the two-dimensional representation of the complex imaginary unit. Indeed, it is well known 
that complex multiplication can be written as matrix multiplication within the group of matrices of the form 

_ \ \ f cos 9 sin9 \ _ 1' Kez Im ^ \ , . 

^ ~ y-sine cose ) ^ y-liRz Rez ) ' ^^^^ 
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where z = |z|exp(i0) is a complex number and the imaginary unit has \z\ = 1 and 6 = 7r/2. Division algebras 
("matrix multiplication algebras") without zero divisors only exist in 1, 2, 4 and 8 dimensions [50, 51]. 

We have discussed the description of a tachyonic wave packet using a pseudo-Hermitian Hamiltonian whose plane- 
wave solutions fulfill the tachyonic dispersion relation. The Hamiltonians iJs = a • p + /3 7^ m, or alternatively 
H' = a ■ p + i 13 m, have pseudo-Hermitian as well as quasi-pseudo-Hermitian properties, as discussed in Sec. Ill A. 
The description leads to a unitary time evolution within wave packets that are composed exclusively of positive- 
energy and negative-energy waves. The scalar product employed in the calculation of the unitary time evolution is 
the ordinary scalar product, and it is positive-definite for positive-energy and negative-energy wave packets. Positive- 
energy solutions are commonly interpreted as those characterizing particles, and negative-energy solution characterize 
antiparticles. It has recently been shown [1] that for neutrinos, the field-theoretical generalization of the derivation 
presented here leads to a natural suppression of right-handed neutrinos and left-handed anti-neutrinos, due to an 
indefinite Hilbert-space norm of the "wrong-helicity" states which is induced by helicity-dependent anticommutators, 
as explained in Ref. [1]. On the level of quantum mechanics, the construction of a conserved (under time evolution) 
scalar product applicable to both positive-energy as well as negative-energy wave packets is left as an open problem. 
However, we stress that due to the ie prescription as employed in Ref. [1] in the tachyonic propagator denominator, 
positive-energy and negative-energy solutions propagate independently. 

Inconclusive indications for a possibly superluminal nature of the neutrino come from low-energy beta decay ex- 
periments [52-58], from an early arrival of neutrinos from the 1987 supernova [59], from the MINOS laboratory-based 
time-of-fiight experiment (as reported in Ref. [60]), as well as from an earlier FERMILAB experiment [61]. The 
OPERA experiment [62] is currently being reexamined. While none of these experiments provides conclusive evidence 
for the superluminal nature of the neutrino, it is interesting to observe that the experimentally determined best esti- 
mate of the neutrino propagation velocity has been measured as superluminal in Refs. [60-62], and the best estimates 
for the neutrino mass squares have been determined to be negative in Refs. [52-58], with a tendency toward smaller 
absolute values for the mass squares at lower energies (for a summary see also [63]). Conceivable energy-dependent 
neutrino mass running is discussed in Refs. [64, 65]. 

Appendix A: Unitary Transformations of the Superluminal Hamiltonian 

In Refs. [30, 37], the Hamiltonian 

H' = d-p + if3m, (Al) 

is investigated, where i is the imaginary unit and m is a real mass term. This Hamiltonian is related to the tachyonic 
Dirac Hamiltonian as follows. We consider the unitary transformation 



U = 



1 



7^ 



( ' 


















1 

















1 






V ^71 





1 

V2 


/ 



(A2) 



which has the property 

[//3 7^C/-i =i/3, a' = Ua'U~\ {a\ a^} ^ {a\ a^} = 2S'^ . (A3) 

Now, 

H" = U H^U-^ (a^p^ + p-f^m) U'^ =a'p' + ipm. (A4) 

Because the 3* fulfill the same algebraic relations as the a*, the Hamiltonian H" is equivalent to the Hamiltonian 
(5* -)■ a') 

H' = a' p' + iPm. (A5) 

Using a further unitary transformation, we may invert the sign of the imaginary mass term, 

H'" = -f^ H' (j^)-^ ^ a'p' ~il3m, (A6) 

where = /3 is a unitary operator. Because of the imitary relation of H'" , H", H' , and H5, it is not surprising 
that the field theories induced by the Hamiltonians are indeed equivalent, as analyzed in detail in Ref. [30] (with a 
special emphasis on i/s and H'). 



15 



Appendix B: Quantum Mechanics and Scalar Products 

We attempt to investigate the connection to quantum mechanical Hamihonians as described in Refs. [10-14]. For 
dcfiniteness, we follow the conventions of Rcf. [43], which implies that the operators V and T are restricted to the 
domain of ordinary (nonrelativistic) quantum mechanical operators, and from now one, C is redefined to he a pseudo 
charge-conjugation in accordance with Ref. [43]. Hamiltonians of the following structure are being investigated, 

H =p^ + x^ {ixY ^—^+x^ (ixY , (Bl) 

where p = —id/dx is the momentum operator. In Refs. [10-14, 43], the parity operator is defined to be equivalent to 
the transformation x — > —x, p — >■ — p, while the T operator describes the time inversion, which is equivalent to a; — > x, 
p — > —p, and i — > — i. We note that the transformation i — > — i is equivalent to the calculation of the Hermitian adjoint. 
For the class of Hamiltonians given by Eq. (Bl), which we consider exclusively from now on, the T transformation is 
equivalent to the calculation of the Hermitian adjoint, and we have 

[VT,H]=Q, H = VT HT^^V^^ =V H+V^^ . (B2) 

In the case of unbroken VT symmetry, the eigenfunctions of the Hamiltonian are also eigenfunctions of the VT 
operator [43], 

H^n{x)^En(t>n{x), V T 4>n{x) = 4>n{x) ■ (B3) 

Using the Hermiticity of the V operator, we can formulate the condition 

^ j Ax [VTfix)] g{x) = J AxTf{x)Vg{x) = j dx r{x)Vg{x) , (B4) 

which is equivalent to our V scalar product discussed in Sec. HI D. Now, for the case of a discrete spectrum, one uses 
Eq. (B3) and obtains [43] 

{4>nArn)= j Ax (j)*^{-x) (l)m{x) = j Ax (j)n{x) (j)m{x) = {-lY 5rim ■ (B5) 

A scalar product is positive definite if its matrix representation is positive definite in a complete basis set of functions. 
Let us choose the basis set to be the set of the eigenstates 0,1 (a;). Then, if we can find an operator C which has the 
following properties, 

C0„(x) = (-l)"0„(x), {f\g) = j Ax[CVTf{x)]g{x), (<^„[ 0™) = , (B6) 

the scalar product ( •[ •) is positive definite. This defines the C operator whose calculation is discussed in Refs. [43, 44]. 

The above formalism is somewhat tied to systems with a discrete spectrum of bound states. For the continuum 
eigenstates of the tachyonic Dirac Hamiltonian, the norm in the V scalar product vanishes because for plane-wave 
states of wave vector fc, it is proportional to the expression 

/ A^x [exp(i k ■ X)] * V exp(i A? • f ) = / d3. cxp(-i • x) exp(.i fc • x) 

= (27r)3,53(2fc) = 0, fcT^O. (B7) 

The physically interesting states have a nonvanishing wave vector and thus, zero V norm. A possible choice for the C 
operator, which is physically consistent within the manifold of positive-energy and negative-energy states, therefore 
is C = P, as discussed in Sec. HID. By "physically consistent" we mean that it leads to a conserved scalar product 
within the time evolution for the positive- and negative-energy states. 
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